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ABSTRACT 



The temperatures of interstellar dust grains are analyzed using stochastic simula- 
tions, taking into account the grain composition and size and the discreteness of the 
photon flux. Grains of radii smaller than about 0.02 are found to exhibit large 
temperature fluctuations with narrow spikes following the absorption of UV photons. 
The temperatures of such grains may rise by a few x 10 K for very short times, but 
they do not rise above 80 K even for irradiation field intensities of photon dominated 
regions. The distribution of grain temperatures is calculated for a broad range of grain 
sizes and for different intensities of the interstellar radiation field, relevant to diffuse 
clouds and to photon dominated regions. The dependence of the average grain tem- 
perature on its size is shown for different irradiation intensities. It is found that the 
average temperatures of grains with radii smaller than about 0.02 fim are reduced due 
to the fluctuations. The average temperatures of grains of radii larger than about 0.35 
/xui are also slightly reduced due to their more efficient emission of infrared radiation, 
particularly when exposed to high irradiation intensities. The average temperatures 
(T) of silicate and carbonaceous grains are found to depend on the radiation field in- 
tensity xmmp according to (T) ~ Xmmp 7 , where the exponent 7 depends on the grain 
size and composition. This fitting formula is expected to be useful in simulations of 
interstellar processes, and can replace commonly used approximations which do not 
account for the grain temperature fluctuations and for the detailed properties of inter- 
stellar dust particles. The implications of the results on molecular hydrogen formation 
are also discussed. It is concluded that grain-temperature fluctuations tend to reduce 
the formation rate of molecular hydrogen, and cannot account for the observations of 
H2 in photon-dominated regions, even in the presence of chemisorption sites. 
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1 INTRODUCTION 

Interstellar dust grains, which consist of silicates and 
carbonaceous materials, account for about 1% of the 
mass of interstellar clouds. The grains play an impor- 
tant role in processes in the interstellar medium (ISM), 
such as the absorption of UV radiation, emission of 
infrared radiation and heating of the surrounding gas 
through the photoelectric effect, which leads to the emis- 
sion of energetic electrons. Catalytic processes on grain 
surfa ces give rise to the fo r mation of molecular hydro- 
gen dGould fc Salpeterlll963l; iHollenbach fc Salpeter] Il97ll ; 



Hollenbach, Werner & Salpeter 1971) and other molecules 
(|Hasegawa. Herbst fc Leung [1992). Dust grains also take 
part in the dynamics of interstellar clouds through their role 
in coupling the magnetic field to the gas in regions of low 
fractional ionization and t hrough their transferring of radia- 
tion pressure into the gas (Drainc 2003). Many of these pro- 



cesses depend directly or indirectly on various properties of 
the grains, such as their temperature and charge. These, in 
turn, depend on both the macro-conditions, namely the ra- 
diation flux and the gas environment, and the microphysics 
of the grains namely their composition and structure. Thus, 
the evolution of interstellar clouds involves a complex net- 
work of physical and chemical processes and feedback mech- 
anisms. Large scale simulations are required in order to com- 
bine all the known processes and to draw the macro-physical 
picture l|van Hooj|2004l : iLe Petit et alj|2006j ). An important 
ingredient of interstellar processes is the temperature of dust 
grains and its dependence on the interstellar conditions. Pro- 
cesses such as infrared emission from the grains and gas- 
grain heating directly affect the grain temperatures. Other 
processes such as the formation of H2 and other molecules on 
grain surfaces are highly sensitive to the grain temperatures. 

In this paper we calculate the temperatures of interstel- 
lar dust grains and their dependence on the radiation field 
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and on the compositions and sizes of the grains. To this end 
we perform simulations of stochastic heating and radiative 
cooling, taking into account the discrete nature of the UV 
photons. We find that grains of radii smaller than 0.2 /xm 
exhibit large temperature fluctuations. Such small grains ex- 
hibit skewed temperature distributions, where most of the 
time they are colder than large grains exposed to similar 
radiation intensity, but occasionally they experience a sharp 
temperature spike. 

The paper is organized as follows. In section[2] we briefly 
review the physical properties of the radiation field in differ- 
ent interstellar environments. Previous studies of interstellar 
grain temperatures are reviewed in section [3] The methods 
used in our stochastic simulations are described in section 
13] The results are presented in section [S] followed by a dis- 
cussion in section [6] and a summary in section [7] 
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2 THE INTERSTELLAR RADIATION FIELD 

The interstellar radiation field is defined as the average ra- 
diation field between stars. There is a widely accepted ap- 
proximation ia£jhet^^c^_S2^xmn_^ the ra diation field, 
proposed bv lMathis. Mezger fc Panagial (|l983l ). Within this 
approximation, the energy density u\ (eV cm -3 /im _1 ) of 
the ambient radiation vs. the wavelength A is given by 
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where xmmp is the enhancement of the starlight com- 
ponent relative to the solar nei ghborhood (|Habing| Il968l ; 
iMathis. Mezger fc Panagial Il983t) and c (cm s" 1 ) is the 
speed of light. The term u x accounts for the UV com- 
po nent of the solar neighborhood. It s explicit form is given 
in iMezger. Mathis fc Panagial (|l983f ). The parameters are 
W 2 = HT 14 , W 3 = HT 13 , W A = 4 x 10" 13 , T 2 = 7500K, 
T 3 = 4000K and T 4 = 3000K, while T = 2.7K is the tem- 
perature of the cosmic background radiation. The function 
B\ (T) is the Planck distribution for black-body radiation at 
temperature T (K) given by 
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where h (eV s) is the Planck constant. The photon flux 
/(A) = \u\/h (photons cm" 2 s" 1 /im" 1 ), obtained from 
Eq. for xmmp = 1, is shown in Fig. [TJ for the relevant 
range of wavelengths 0.0912 < A < 975 /xm (or photon en- 
ergies in the range 1.27 x 10" 3 < E < 13.6 eV). Photons 
that are more energetic than 13.6 eV are suppressed inside 
the cloud because they ionize hydrogen atoms on the edge 
of the cloud and cannot penetrate further. 



3 DUST-GRAIN TEMPERATURES 

Information on the properties of interstellar dust grains 
is obtained from spectroscopic observations of the extinc- 
tion, scattering and emission of electromagnetic radiation by 
these grains. Such observations show that the grains consist 
of silicate and carbonaceous materials. Silicate grains consist 
of (Mg^Fei-^SiCU, where < x < 1. Carbonaceous grains 
include crystalline graphitic carbon, amorphous carbon and 



Figure 1. The interstellar radiation field /(A), expressed in terms 
of the number of photons per unit wavelength (cm) , per cm 2 , per 
second, for xmmp = 1- 



polycyclic aromatic hydrocarbons (PAHs). In cold molecu- 
lar clouds the grains are coated by ice mantles. However, 
here we focus on bare silicate and carbon grains, relevant to 
diffuse clouds and photon-dominated regions. 

The size distribution of interstellar dust grains fol- 
lows an approximate power law of the form dn(r)/dr oc 
r _A , where r is the grain radius, n(r) (cm -3 ) is the den- 
sity of grains smaller than r and A = 2.5. Such dis- 
tribution, between suitable cutoffs r m i n < r < r max , 
was shown to reproduce the ob served extinction curve 
l|Mathis. Rumpl fc NordsieckHl977l ). Later work by Draine 
& Lee (1984) provided further evidence for the model. How- 
ever a more recent study by Weingartner & Draine (2001) 
has shown that a range of different size distributions ex- 
ist compatible with observations of different regions in the 
Milky Way, the LMC and the SMC galaxies. 

The temperatures of interstellar dust grains depend on 
the heating and cooling processes due to the interaction of 
these grains with the radiation fields and the surrounding 
gas. These processes include the absorption and emission 
of electromagnetic radiation, gas-grain collisions, heating by 
cosmic rays and surface reactions. In most interstellar envi- 
ronments, the grain heating and cooling processes are dom- 
inated by absorption and emission of radiation. Gas-grain 
collisions strongly affect grain temperatures only i n dense 
regio ns where the rate of such collisions is high l|Spitzer] 
Il978h . Chemical reactions affect the grain temperatures in 
very cold dense regions in the cores of molecular clouds. 

Dust grain temperatures have been studied for many 
years. DeHulst (1949) estimated the temperature of di- 
electric particles to be ~ 15K. Detailed calculations 
of the temperatures of grains of different co mpositions 

in different radiation fields were done by iGreenberg 

dl97ll). Later, a graphite-si li cate model was proposed by 
IMathis. Rumpl fc Nordsieckl l| 19771 ). 

The t emperature fluctuati ons of s mall grains w ere an- 
alyzed by iDraine fc Led |l984h and by iDwekl (|l986l ). using 
statistical methods. These focused mainly on the infrared 
emission from such grains. Li & Draine (2001) evaluated the 
temperatures of dust grains as a function of their sizes and 
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the ambient radiation field, using the balance between the 
absorbed and emitted radiation: 

/ Qnbs(r, \)cu\d\ = / iTrQ &h3 (r,\)B x (T)d\. (3) 
Jo Jo 

Here, vTr 2 Q a bs(r, A) is the absorption cross-section for a pho- 
ton of wavelength A (cm) by a grain of radius r where 
Qabs(r, A) is the absorption coefficient. The energy density 
vs. wavelength, u\ (eV cm -3 fim" 1 ), is given by Eq. ([1]). 
Thus, cu\ (eV cm -2 s" 1 /im" 1 ) is the energy flux density 
at wavelength A. Eq. Q describes the balance between ab- 
sorption (left-hand side) and emission (right-hand side) for 
large grains that maintain a nearly constant temperature T 
(K). While most of the energy absorbed by grains is in the 
UV range, the emitted radiation is in the infrared range. 
From Eq. ([3]) one can extract the grain temperature T vs. 
the grain size r and the radiation field intensity, given by 
Xmmp- Note that the temperatures of such large grains do 
not depend on the heat capacity of the grain. 
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Figure 2. The cross-section Q a b s (A) for photon absorption onto 
a graphite grain, as a function of the wavelength. 



4 THE STOCHASTIC SIMULATIONS 

Here we present stochastic simulations of grain heating and 
cooling, taking into account the grain composition and size 
and the discrete nature of the UV photons. In the simula- 
tions we follow the temporal evolution of the temperature 
of a single grain. The spectrum of absorbed photons is di- 
vided into two domains: (a) the absorption of an energetic 
UV photon is taken as a discrete, stochastic process, (b) The 
contribution of the less energetic photons is considered as a 
continuous flux of energy. The emission of infrared radiation 
by the grains, is also described as a continuous process. 

The time-dependent temperature of a single grain is 
evaluated as follows. Starting from a given initial tempera- 
ture, the grain temperature is updated taking into account 
both the discrete and the continuous processes. The contin- 
uous variation of T(t) is calculated by integrating the ab- 
sorbed and emitted radiation, taking into account the heat 
capacity of the grain. To evaluate the stochastic contribu- 
tion we specify the times in which UV photons are absorbed 
by the grain. These times are drawn from a suitable Pois- 
son distribution, which accounts for the flux of energetic 
photons and the wavelength-dependent cross section of the 
grain. For each absorption event, the wavelength of the pho- 
ton is drawn from a suitable distrituion obtained from u\ 
and Qabs- The grain temperature is then updated according 
to the energy of the photon and the temperature-dependent 
heat capacity of the grain. 



4.1 Physical Properties of Dust Grains 

For simplicity, we assume that the grains are spherical and 
denote the grain radius by r. Two grain compositions are 
studied: carbon (or graphite) grains with mass density of 
2.16 (gram cm -3 ) and silicate grains with mass density of 3.5 
(gram cm -3 ). In the analysis we use the physical properties 
of such spherical, single-component grains given by Laor & 
Draine (1993). The photon absorption coefficient Q a b s as 
function of the wavelength A for graphite grains is shown in 
Fig. [2] This absorption coefficient depends on the grain size. 
It is also weakly dependent on the temperature, however this 
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Figure 3. The heat capacity per carbon atom in the graphite 
grain (solid line) and per atom in the silicate grain (dashed line) 
vs. grain temperature. Here C(T) stands for C grap hite(T) or 
C s iii ca te(T), while N is N c or N a accordingly (Draine et al. 2001). 
the Boltzmann constant kg is given in units of eV/K. 



dependence is neglected in our analysis. The values of Q a bs 
used here are for grain temperature of 25 K (Laor & Draine 
1993). 

The heat capacities of silicate and graphite grains are 
shown in Fig. [3] The heat capacity of a graphite grai n at 
temperature T can be expressed by (|Draine fc Lill200ll ) 



C graphite (T) = (N-2)ki 



f ' 2 (lis) + 2/2 (2504) 



(4) 



where TV is the number of carbon atoms in the grain and 
f' n {x) = df n {x)/dx, where 
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Similarly, the heat capacity of a silicate grain is given by 
i|Draine fe Lill200lD 



C si iic atc (T) = (N- 2)ki 



2/2 (50o) +/3 ( 
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1500 



(6) 



where iV is the number of atoms in the grain. 
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4.2 The Photon Flux 

In order to simulate the time dependent temperature of a 
grain exposed to a given radiation field, we first divide the 
relevant spectrum of the radiation field into two domains: 
the hard, energetic photons, taken discretely, and the soft 
photons, taken as a continuous flux. For the discrete part 
we evaluate the average time t between absorptions of hard 
photons. The actual absorption times are then obtained, us- 
ing a poisson distribution, with an average r, for the times 
between absorption events. The wavelength of each absorbed 
photon is drawn from a suitable distribution, which accounts 
for the interstellar radiation field and the absorption coeffi- 
cient Qabs- 

The threshold A c (r) between the domains of hard and 
soft photons depends on the grain radius r. The criterion for 
hard photons is that these photons are sufficiently energetic 
to cause a significant change in the grain temperature upon 
absorption. In practice, A c (r) is determined such that pho- 
tons with A < A c raise the temperature of a grain of radius 
r from To to T > T + AT. More specifically we use T = 10 
K and AT = 0.005 K. These values provide a good descrip- 
tion of the stochastic features while maintaining a reasonable 
running time of the simulation even for large grains. 



4.3 Stochastic Heating of Grains 

The stochastic heating of a grain of radius r exposed to radi- 
ation field u\ is simulated as follows. First, we evaluate the 
average time, r, between succes sive absorptions of energet ic 
photons. This time is given by (|Draine fc Andersonl ll985') 



X e (r) 



/(A)Q abs (r,A)dA 
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where A m in = 0.0912 (fim) is the lower cutoff of the spec- 
trum, which corresponds to ionization threshold of atomic 
hydrogen. The times, At, between photon absorption events 
are drawn from the distribution P(At) = exp(— At/r)/r. 
The next step is to obtain the wavelength (and thus the 
energy) of the absorbed photon. This wavelength is drawn 
from the distribution 



P(A) = 



i(A)Q abs (r,A) 
f A Ac(r) 7(A)Q abs (r, A)dA 



(8) 



The temperature T of a grain immediately after the 
absorption of a UV photon of energy -Bphoton = hc/X is 
calculated using the relation 
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where C(T') is the heat capacity of the grain and To is the 
instantaneous grain temperature before the absorption of 
the photon. To avoid an accumulation of numerical errors, 
the temperature T is calculated analytically. To this end, 
the heat capacity C(T) is fitted by a polynomial function, 
which is then integrated, giving rise to an expression of T 
in terms of To and .Ephoton- 



a 4o 



(a) 



t— 10.000 20.000 30.000 40.000 5,0000 6.0000 7.0000 80.000 90,000 100,000 



(b) 



t— 10.000 20,000 30.000 40,000 50,000 60.000 70.000 80.000 90,000 100,000 

--80 
| 40 



t— 10.000 20,000 30,000 40,000 50,000 60,000 70,000 80.000 90,000 100.000 
Time (s) 



Figure 4. The temperature variations of a carbon grain during 
one day in which it is exposed radiation intensity of Xmmp = 1; 
where the grain radii are: (a) 0.05; (b) 0.016; and (c) 0.005 fim. 



4.4 Continuous Heating and Radiative Cooling of 
Grains 

The time dependence of the grain temperature between ab- 
sorptions of UV photons is determined by the balance be- 
tween the absorption of soft photons and the emitted radi- 
ation. The time derivative of the grain temperature is thus 
given by 



dT 
dt 

where 
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(10) 



(11) 



is the radiation absorbed by a grain of radius r within 
the frequency range of soft photons. To obtain the time- 
dependent grain temperature T(t), Eq. (|10fl is integrated 
using a standard Runge Kutta stepper during the period 
between succesive absorption events of hard photons. 



5 RESULTS 

5.1 Grain Temperature vs. Time 

In Fig. [4] we present the grain temperature vs. time for car- 
bon grains with radii of (a) 0.05, (b) 0.016 and (c) 0.005 fim 
for radiation intensity of Xmmp = 1. Clearly, for the largest 
grain size the temperature fluctuations are negligible while 
for the smallest grain they are intense. The temperature of 
the small grains is low most of the time (10-15K). However, 
following the absorption of a UV photon, the grain tem- 
perature may rise up to nearly 80K followed by fast cool- 
ing, giving rise to narrow temperature spikes. This is due 
to the fact that the emitted radiation intensity sharply in- 
creases as a function of the grain temperature. As a result, 
the fluctuations tend to reduce the average temperature of 
the small grains. In Fig. [5] we present the time dependent 
grain temperatures with a higher temporal resolution. The 
dependence of the rate of UV photon absorptions on the 
grain size is clearly observed. 
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Figure 5. The temperature variation of the three carbon grains 
shown in Fig. [4] during 1000 seconds in which they are exposed 
to radiation intensity of Xmmp = 1- 
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Figure 7. The average temperature and fluctuation range vs. 
grain radius for carbon grains exposed to different irradiation field 
intensities. 
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Figure 6. The distributions of grain tempeatures for grains of 
radii of 0.005 (solid line), 0.02 (dashed line) and 0.08 (dashed- 
dotted line) fim: (a) carbon grain exposed to radiation of Xmmp = 
1; (b) carbon grain exposed to Xmmp = 100; (c) silicate grain 
exposed to xmmp = 1; and (d) silicate grain exposed to xmmp = 
100. 



5.2 Grain- Temperature Distributions 

Due to the wide variations in the temperatures of small 
grains, one should calculate not only the average temper- 
ature (T) but the entire probability density function f(T). 
This density function is obtained directly from T(t) by cal- 
culating what fraction of time the grain has spent in each 
range of temperatures. In Fig. [5] we present the distribution 
of grain temperatures for carbon grains [Figs.[6ja) and[6jb)] 
and for silicate grains [Figs. |6jc) and[6jd)]. The results are 
shown for three grain radii, 5 (solid line), 20 (dashed line) 
and 80 nm (dashed-dotted line) and two radiation field in- 
tensities: Xmmp = 1 [Figs. [6ja) and |6jc)] and Xmmp = 100 
[Figs. ®h) andHJJd)]. The large grains exhibit very narrow 
temperature distributions. This is due to the fact that the 
large grains absorb a high flux of UV photons and that the 
effect of each photon on their temperature is small. As the 
radiation field increases the peak shifts to higher temper- 
atures. The small grains exhibit broad temperature distri- 
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Figure 8. The average temperature and fluctuation range vs. 
grain radius for silicate grains exposed to different irradiation field 
intensities. 



butions. These distributions are skewed, where most of the 
weight is concentrated at the low temperature side with only 
a narrow tail at the high temperature side. The typical tem- 
perature of the small grains is lower by about 10 K compared 
to the large grains. 



5.3 Average Grain Temperatures 

In Figs. [7] and [5] we present the average grain tempera- 
tures (T) (solid lines) vs. grain radius for carbon and silicate 
grains, respectively, for radiation intensities of xmmp = 1, 
10, 100 and 1000. The gray areas represent the range of vari- 
ation of the grain temperature. The upper (lower) bounds, 
T+(r) (T_(r)) of these areas represent the upper (lower) 
semi-standard deviations with respect to (T), computed 
from the probability density function f(T), taking into ac- 
count only the range T > (T) (T < (T)). More specifically, 



{T) + 
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and 
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Table 1. Fitting parameters for the average grain temperature 
vs. radiation field for carbon grains 



Table 2. Fitting parameters for the average grain temperature 
vs. radiation field for silicate grains 



Grain Radius (fJ.m) K 



Grain Radius (|im) K 



0.005 14.5 0.23 

0.05 21.0 0.18 

0.5 16.4 0.16 



0.005 12.6 0.21 

0.05 16.8 0.16 

0.5 15.2 0.16 



T_(r) = (T> 



So \T-{T) ff{T) 
So T) f(T) 



1/2 



(13) 



It is found that for all radiation fields, grains of radii larger 
than r c — 0.02 /an do not exhibit significant temperature 
fluctuations. For grains of radii r < r c , the temperature 
fluctuations are enhanced as d is reduced. The threshold 
value of the diameter below which fluctuations are signif- 
icant is clearly independent of the radiation field. This is 
due to the fact that the fluctuation level is determined by 
the temperature rise that is caused by a single UV photon, 
which is reduced as the grain mass is increased. For small 
grains, the temperature spikes caused by UV photons are 
very narrow and exhibit little overlap even for high radia- 
tion fields. For high radiation fields the average temperature, 
(T) (solid line) is at about the middle of the gray area. For 
small grains, the average temperature is near the lower edge 
of the gray area. This is due to the fact that the small grains 
spend most of the time at low temperatures. 

The results for the temperature distributions (Fig. 6) 
and for the average grain temperatures (Figs. 7 and 8), for 
Xmmp = 1, are in agreement with Li & Draine (2001) and 
with Cuppen, Morata & Herbst (2006) in the correspond- 
ing range of grain sizes, up to variations of the order of 
IK. The results presented here extend the analysis of grain 
temperatures to higher radiation fields. These results can 
thus be used in modeling of both diffuse clouds and photon- 
dominated regions. 

The highest average temperatures are obtained for 
grains of intermediate size, with radii in the range of 20-35 
nm. Smaller and larger grains exhibit lower average tem- 
peratures. This can be understood as follows. Small grains 
radiate much more efficiently during the short times they 
spend at high temperatures and thus tend to have lower 
average temperatures. Also, very small grains absorb UV 
radiation less efficiently than large grains, in both graphite 
(Fig. 2) and silicate grains. Large grains radiate much more 
efficiently in the relevant range of infrared radiation. The 
average temperatures of grains of different sizes as a func- 
tion of the radiation field intensity xmmp are shown in Figs. 
[9] and [10] . These curves are fitted by the function 



(T) = A"xmmp" 



(14) 



where K and 7 depend on the grain size. The values of K 
and 7 for grains of three different sizes, which consist of 
silicates and carbon are shown in Tables 1 and 2. 

To interpret the resulting values of 7, consider first a 
macroscopic object which absorbs and emits radiation like 
a black body. For such object Q a b s = 1. Its temperature T 
does not exhibit fluctuations. Since the intensity of the ra- 
diation emitted from such object is proportional to T 4 , one 




Figure 9. The average temperature vs. interstellar radiation 
field Xmmp f° r carbon grains with radii of 0.005 (solid line), 0.05 
(dashed line) and 0.5 (dashed-dotted line) fim. 



obtains that 7 = 1/4. Unlike a black body, for large interstel- 
lar grains, in the infrared regime, Q a bs ~ A~ 2 . According to 
Wien's law, for an object at temperature T, the highest radi- 
ation flux is emitted at wavelength A max (T) = w/T, where 
w — 0.29 (cm K) is the Wien's constant. As a result, the 
emitted radiation is proportional to Q a b s (A ma x(T'))T 4 ~ T 6 , 
while the incoming radiation is proportional to xmmp- Thus, 
T ~ xmmp 1 ' 6 , namely 7 = 1/6. This analysis agrees very 
well with the results shown in Tables 1 and 2 for silicate and 
carbon grains of radii r — 0.5 /im. For very small grains, 
stochastic fluctuations reduce the average temperature in 
the limit of low radiation intensity. As a result, larger values 
of 7 are obtained. 




Figure 10. The average temperature vs. interstellar radiation 
field Xmmp f° r silicate grains with radii of 0.005 (solid line), 0.05 
(dashed line) and 0.5 (dashed-dotted line) fim. 
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6 DISCUSSION 

The results presented above are expected to be useful in sim- 
ulations of processes taking place in interstellar clouds. In 
particular, the average grain temperatures vs. grain size and 
radiation intensity can be used as input data for such simu- 
lations. Also, the grain-temperature distribution, combined 
with the grain size distribution can be used to evaluate the 
infrared radiation spectrum emitted from dust grains in an 
interstellar cloud. 

The commonly used approach to the evaluation of grain 
temperatures are based on the balance between the absorbed 
and emitted radiation, given by Eq. (3). Using Eq. (3) as an 
implicit equation, the grain temperature can be extracted 
numerically. In this case one can take into account the de- 
pendence of the absorption and emission properties on the 
grain size, through the function Q a b s - 

In this approach the stochastic heating which leads to 
temperature fluctuations of small grains is ignored. Note 
that due to the broad distribution of grain sizes, most of 
the surface area of interstellar dust is in small grains. Thus 
small grains dominate the radiation emitted from grains. 
In order to obtain an explicit formula for the grain tem- 
perature, the function Qabs is often replaced by a crude 
approximation, in which it is constant in the UV range, 
and decays as 1/A in the infrared ran ge [see Eq. (2) in 
iHollenbach. Takahashi fc Tielensl i|l99ll )]. This approxima- 
tion ingnores the specific features of the silicate and carbon 
materials as well as the dependence of Q a bs on the grain 
size. For large grains, the assumption of Q a bs = 1 in the 
UV range is valid. However, in Fig. 2 it appears that in 
the range of wavelengths in which most of the radiation 
from grains takes place (A > 30 fim) the functional de- 
pendence is Q a bs ~ 1/A 2 rather than 1 / A as assumed in 
IHollenbach. Takahashi fc Tielensl (|l99ll) ]. This feature di- 
rectly affects the dependence of the average grain tempera- 
ture on the radiation field intensity (see Sec. 5.3 above). For 
small grains, the assumption of constant Q a b s is not valid. 
More specifically, for small grains Qabs strongly depends on 
A and on the grain size (Fig. 2). 

Our analysis takes into account both the specific fea- 
tures of the silicate and carbon materials, including the de- 
tailed dependence of the absorption and emission on the 
wavelength and on the grain size. In addition, the effect 
of stochastic heating is taken into account. This approach 
provides a more precise evaluation of interestellar dust tem- 
peratures under a broad range of conditions. It shows the 
dependence of the average grain temperature on both the 
grain composition and size and on the radiation field inten- 
sity. 

Experimental studies show that molecular hydrogen for- 
mation on amorphous silicate and carbon grains is efficient 
only in a narr ow temperature windo w in the ra n ge be - 
tween 10-20 K iPirronello et all i|l997h ; iKatz et al.l l)l999h : 
iPerets et all (|2007l ). This range coincides with the typical 
grain temperatures in diffuse clouds. Therefore, such grains 
are efficient catalysts for H2 formation in diffuse clouds, as 
long as they are sufficiently large. The efficiency of molecular 
hydrogen formation on small grains is reduced under condi- 
tions in which the nu mber of adsorbed H atoms on the grain 
is of order 1 or lower (iTielens fc Hagenlll982l ; ICharnlev et al.l 
Il997l ; I Caselli et al.lll998h . Under these conditions, the rate 



equations fail and stochastic methods based on the master 
equation are needed in order to evaluate the reactio n rates 
l|Biham et al.ll200ll ; iGreen et al1l200ll ; ICharnlevll200ll) . Since 
the grain size distribution is dominated by small grains, their 
effect on the rate of molecular hydrogen formation is impor- 
tant. 

The experimental results cannot explain the high abun- 
dance of H2 in photon-dominat ed regions wh e re th e grain 
temperatures may reach 50 K lHabart et al.l l|2004h . This 
led to the possibility that temperature fluctuations of small 
grains may enhance the formation rate of molecular hydro- 
gen and possibly resolve the discrepancy between the labo- 
ratory experiments and astrophysical observations. 

The effect of grain temperature fluctuations on 
the formation rate of molecular hydrogen on small 
dust grains in diffuse cloud s was studied recently 
|Cuppen. Morata fc Hcrbst 2006). It was found that for ra- 
diation intensity characterized by xmmp = 1, molecular hy- 
drogen formation is efficient only on very rough grains which 
exhibit higher energy barr iers for H desorpt io n than those 
measu red experimentally Katz et al. I (|l999h : IPerets et al.l 
(2007). However, deeper in the cloud, where the radiation 
intensity is reduced by extinction, molecular hydrogen for- 
mation is efficient also using the experimentally measured 
grain parameters. 

To understand the effect of stochastic heating on hydro- 
gen recombination, consider a small grain exposed to UV ra- 
diation. During the time between the temperature spikes the 
grain is colder than a large grain exposed to a similar radi- 
ation field. If the time between spikes is long enough for the 
grain to adsorb several H atoms, these atoms may quickly re- 
combine during the next spike. In this case, stochastic heat- 
ing may enhance H2 formation on small grains. However, 
even for xmmp = 1 the rate of UV photon absorption on a 
grain is much higher than the adsorption rate of H atoms. 
As a result, small grains rarely accumulate more than one H 
atom between temperature spikes. Thus, stochastic heating 
alone does not enhance the rate of H2 formation. 

It was suggested that chemisorption sites on grain sur- 
faces may also en hance the formation of H2 in photon- 
dominated regions l|Cazaux fc Tielens|[2002l . l2004r ). However, 
experiments on H2 formation on graphite indicate that in 
order to enter the chemisorption sites the H atoms may 
need to pass through an energy barrier (|Zecho et al.l 120021 ; 
IZecho. Guttler fc Kuppersll2004r i. In addition, molecular hy- 
drogen formation from chemisorption sites is markedly dif- 
ferent than the Langmuir-Hinshelwood mechanism observed 
in physisorption sites. In particular, H2 molecules are formed 
only upon desorption of the H atoms from the surface, 
which occurs only at high surface temperatures of T > 
400 K JZecho et al.l |2002| ; IZecho. Guttler fc Kuppersl |2004| ; 
IPerets fc BiharrfeOOrj T These temperatures are much higher 
than the maximal grain temperatures of stochastically- 
heated small grains in photon-dominated regions (Fig. 6). 
Therefore, it seems that even the combination of stochastic 
heating and chemisorption does not account for the forma- 
tion of H2 in photon-dominated regions. 
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7 SUMMARY 

The temperatures of interstellar dust grains and their depen- 
dence on the grain size and composition and on the radiation 
field intensity have been studied using stochastic simula- 
tions. It was found that grains of radii smaller than 0.02 /jm 
exhibit temperature fluctuations which become more intense 
as the grain size decreases. The temperature distribution of 
small grains is a skewed function, dominated by the low 
temperature side with a narrow tail in the high temperature 
side. The fluctuations give rise to a reduction in the aver- 
age grain temperature. The average temperatures of grains 
of radii larger than 0.035 /im are also slightly reduced due 
to their more efficient emission of infrared radiation. The 
average grain temperatures vs. radiation intensity were fit- 
ted according to (T) = Kxmmp 1 , where K and 7 depend 
on the grain size and composition. These results are more 
accurate and general than the commonly used formulae for 
interstellar dust temperatures, which do not take into ac- 
count the detailed features of the absorption ane emission 
and the effects of stochastic heating. 
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